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We study extensions of the minimal supersymmetric standard model (MSSM) with new degrees 
of freedom that couple sizably to the MSSM Higgs sector and lie in the TeV range. After integrating 
out the physics at the TeV scale, the resulting Higgs spectrum can significantly differ from typical 
supersymmetric scenarios, thereby providing a window beyond the MSSM (BMSSM). Taking into 
account current LEP and Tevatron constraints, we perform an in-depth analysis of the Higgs collider 
phenomenology and explore distinctive characteristics of our scenario with respect to both the 
standard model and the MSSM. We propose benchmark scenarios to illustrate specific features of 
BMSSM Higgs searches at the Tevatron and the LHC. 



I. INTRODUCTION 

There has been a recent surge of interest in extensions 
of the minimal supersymmetric standard model (MSSM) 
by higher-dimension operators [THU]- These can have 
an important impact on the Higgs sector, alleviating in 
particular the tension present in the MSSM that results 
from the LEP Higgs bounds. Such effective field theory 
(EFT) studies allow a model-independent description of 
a large class of extensions of the MSSM, and permit one 
to quantify the sense in which the Higgs sector can serve 
as a window beyond the MSSM (BMSSM). 

This point of view was clearly put forward in Ref . |3] , 
where it was emphasized that at leading order in 1/M ~ 
where M is the scale of the physics that is integrated out- 
the MSSM is extended by only two parameters. The sur- 
prisingly large effects of such higher-dimension operators 
can be understood from the fact that the MSSM Higgs 
potential is rather restricted at tree-level. The nonrenor- 
malizable operators in the superpotential induce renor- 
malizable (quartic) operators in the Higgs potential that 
are not present in the MSSM limit (at tree- level) , so that 
in spite of the fact that their coefficients are "small" -of 
order fi/M- they correspond to qualitatively new effects. 
In fact, the operators thus induced can easily "destabi- 
lize" the MSSM-like minimum and lead to new minima 
that exist only as a direct result of the higher-dimension 
operators (i.e. the heavy physics). It was emphasized in 
[7] that such minima can be phenomenologically viable, 
can be studied within the EFT framework, and can ex- 
plain the distinct properties induced by the heavy physics 
on the Higgs sector. 

If the BMSSM physics is sufficiently heavy, the lead- 
ing order analysis at order 1/M can suffice. However, 
it is perfectly possible that M is not too far from the 
electroweak (EW) scale, and that nevertheless the heavy 



physics may not be easy to see directly at the LHC, even 
if it is within its kinematic reach (e.g. heavy singlets that 
couple only through the Higgs sector). In such cases, the 
EFT approach is still useful to describe the properties of 
the MSSM Higgs sector. It turns out that the effects of 
order 1/Af^ are more important than naively expected. 
This observation also finds a simple explanation in the 
structure of the MSSM tree-level Higgs potential [8] to- 
gether with the smallness of the MSSM tree-level quartic 
couplings (the root cause for a Higgs state lighter than 
the Z mass in the MSSM at tree-level). The crucial point 
is that the leading order contributions to a subset of the 
quartic Higgs operators, from the heavy physics, first en- 
ter at order 1/A/^. Thus, these end up correcting a coef- 
ficient of order g'^ instead of a 1/M effect, and can give a 
relevant correction even if the expansion parameter is rel- 
atively small. Nevertheless, it is important to appreciate 
that the fact that the first two orders in the expansion 
in 1/M can even result in comparable contributions to 
the Higgs masses, in no way implies a breakdown of the 
EFT. 

In Ref. [H] a detailed study of the consequences for the 
Higgs masses and couplings up to order 1/M^ was given, 
and a selected number of phenonienological observations 
were already made, such as: enhanced gluon fusion pro- 
duction cross sections in a large number of cases, and 
the presence of "exotic" decay modes with more than 
one Higgs boson in the decay chain. In this work we 
analyze the constraints from LEP and the Tevatron on 
the neutral Higgs bosons, as well as the charged Higgs 
bounds from LEP [? ] . We also expand on the associated 
collider phenomenology, emphasizing the type of signals 
that can be expected at both the Tevatron and the LHC. 
We point out that due to the large corrections to the 
Higgs masses (especially to the CP-even Higgs bosons) 
the production and decay patterns can be markedly dif- 
ferent from those in the MSSM. Examples include mod- 



els where both CP even Higgs bosons have significant 
branching fractions into gauge bosons, thus giving rise 
to spectacular signals such as two clearly defined peaks 
in the di-lepton invariant mass distribution. In addition, 
we observe new decay chains that allow for production 
of the "nonstandard" Higgs bosons without large tan /3 
enhancements. It is possible that the full two-Higgs- 
doublet-model (2IIDM) content can be mapped in de- 
tail, thus providing a clear and definite signal for physics 
beyond the standard model, and a rather detailed un- 
derstanding of the mechanism of electroweak symmetry 
breaking (EWSB). If, in addition, relatively light super- 
particle signals are observed, as might be expected in 
these scenarios, a clear case for BMSSM physics could be 
established. Apart from the collider phenomenology in- 
duced indirectly by the heavy physics, higher-dimension 
operators have also been studied in the context of dark 
matter [T2HI4] , cosmology [H] and EW baryogenesis [16l - 
[5D] , and it may be interesting to further explore the con- 
nections with collider physics. 

This paper is organized as follows. In Sec. Ill] we sum- 
marize the most relevant aspects of the models under 
study. In Sec. |III| we discuss the modifications of the 
Higgs spectrum, which are the dominant factor in de- 
termining the Higgs collider phenomenology. In Sec. |IV| 
we discuss in detail the range of signatures uncovered by 
our survey, separating the analysis into the low and large 
tan/3 regimes. We conclude in Sec.|VJ 



II. EXTENDED SUSY HIGGS SECTORS AT A 
GLANCE 

As already mentioned, when considering BMSSM sce- 
narios where the non-MSSM degrees of freedom have 
masses parametrically larger than the weak scale, an 
EFT approach is very useful. The fact that at lead- 
ing order only two parameters are added to those in the 
MSSM makes this a rather economic extension |3], that 
nevertheless can significantly change the MSSM Higgs 
phenomenology. However, the same reason that makes 
these 1/Af suppressed effects rather important also im- 
plies that the next order in the 1/M expansion can be 
phenomenologically relevant, without implying a break- 
down of the EFT [5]. At order 1/M^ there are several 
SUSY-preserving and SUSY- violating operators in the 
Kahler potential, the most important of which, in rela- 
tion to the Higgs phenomenology, were listed in Ref. [8]. 
We refer the reader to this reference for the detailed form 
of such operators and how they affect the expressions for 
the Higgs masses and couplings. Here we restrict our- 
selves to a few general remarks that summarize the most 
relevant features for the present study (full details were 
given in the above reference). 

First, it has to be pointed out that the higher- 
dimension operators to order 1/M^ can be easily gen- 
erated from UV completions that include a combination 
of Higgs singlets, SU{2) Higgs triplets, heavy W primes 



and Z primes. As argued in J5^ the upshot is that the co- 
efficients of the higher-dimension operators, from a low- 
energy point of view, can be chosen in an uncorrelated 
manner. Although the EFT description to order 1/M^ 
introduces a large number of parameters, which makes 
the framework more involved compared to the truncation 
at order 1/Af, one should notice that this same feature 
gives additional handles to infer properties of the heavy 
sector from the properties of the low-energy degrees of 
freedom. In any case, since our goal is to survey the col- 
lider signal possibilities in a model-independent way (in 
a supersymmetric framework), we focus on a low-energy 
study based on the EFT at order 1/A/^, as described 

in m- [? ] 

A random scan over parameter space was performed, 
and a set of points satisfying several constraints was se- 
lected. The set of points in this study satisfy: 

• All the dimensionless coefficients parametrizing the 
higher-dimension operators are taken to be at most 
of order one, i.e. it is assumed that the heavy 
physics at M is weakly coupled. 

• Global minimum: since the scalar potential can 
present several minima, we make sure that the 
vacuum under study is the global one (at least 
within the EFT). We also check that there are no 
charge/color breaking minima, and for simplicity 
we restrict to the CP conserving case (checking 
that the global minimum does not break CP spon- 
taneously). 

• Robustness: there are no accidental cancellations 
that can render (not computed) higher orders in 
the 1/Af expansion more important than expected. 

• "Light" SUSY spectrum: given that generically, 
and unlike in the MSSM, these models satisfy the 
LEP bounds on the Higgs mass at tree-level, there 
is no need for large radiative corrections. Natural- 
ness suggests that in these models the SUSY spec- 
trum would be expected to be light (in the few hun- 
dred GeV range, consistent with direct bounds). 

• Agreement with EW precision constraints, in par- 
ticular in regards to the Peskin-Takeuchi T pa- 
rameter [21]. These arise from three sources: a 
subset of the higher-dimension operators (as gen- 
erated, for instance, by Higgs triplets), the de- 
tails of the MSSM Higgs spectrum, and potential 
custodially- violating mass splittings in the sparticle 
spectrum. We emphasize that mild cancellations 
allow for higher-dimension operator effects that can 
have a non-negligible impact on the Higgs collider 
phenomenology. 

All of the above constraints were described in detail in 
[8j. In addition, we impose the current collider bounds 
from LEP and the Tevatron using the code HiggsBounds 
vl.2.0 mils]. [? ] To this we add the LEP bounds on 



charged Higgs production [53], and the newest combined 
result from the Tevatron in the WW channel [25] , and in 
the inclusive tau search [55] , that are not included in the 
currently available version of this code. We use Higgs- 
Bounds in the "effective coupling" mode, which requires 
effective couplings defined by 
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where (/> = h,H,A is any of the neutral Higgs states, 
T{(j) — >■ X) is the partial width in our model into any 
of the final states X — ss, cc, bb, tt, WW, ZZ, 77 or gg 
(when applicable), and Tsm{4' ~^ -^) is the partial width 
for a SM Higgs of the corresponding mass. Together with 
the total widths in our model (and in the SM), these ef- 
fective couplings encode the information about branching 
fractions into these decay channels in our model. 

We have implemented our tree-level expressions for the 
spectrum and Higgs couplings in HDECAY v3.4 [57] . 
This allows us to compute the Higgs partial decay widths, 
taking into account the QCD radiative corrections, that 
are known to be sizable (for a review, see |28j). In ad- 
dition, we include the radiative corrections derived from 
the 1-loop RG improved effective potential due to super- 
symmetric particles [H], and the SUSY QCD/EW cor- 
rections to the Yukawa couplings [301 EI] • Loop contri- 
butions from the heavy physics that has been integrated 
out are suppressed by both a loop factor and by powers 
of M, hence they are expected to be negligible. 

In all the plots that follow, we have fixed the 
following dimensionful parameters: M = 1 TeV, 
/i = ms = 200 GeV, [? ] and for simplicity, we use 
a common value Msusy = 300 GeV and At = Ab = 
in the MSSM sparticle sector. [? ] The light superparti- 
cle spectrum implies that the loop contributions to the 
Higgs masses are modest, while the loop contributions to 
the Higgs couplings are more important and sensitive to 
the details of this spectrum [301^35]. The above choice 
of Msusy = 300 GeV is simply meant to illustrate the 
possible loop effects arising from relatively light super- 
particles. In particular, one can expect the first two gen- 
eration squarks to be somewhat heavier to satisfy direct 
collider bounds [331 [31] or the sleptons could be some- 
what lighter, without changing our generic conclusions 
regarding the Higgs collider phenomenology. Note also 
that the neutralino/chargino sector depends on parame- 
ters not affecting the Higgs sector directly, and in partic- 
ular that we do not impose constraints from dark matter 
(in this work, we remain agnostic as to the identity of the 
DM candidate, but see I3S]). We have also not imposed 
indirect constraints, such as those arising from b — )• 57, 
Bg — ^ fJ-^ fJ-^ and 3^ — 2, that have the potential to put im- 
portant restrictions, but depend on the flavor structure 
of the soft SUSY-breaking parameters. 

We consider two representative values of tan /3: tan /3 — 
2 and tan/3 = 20. The CP-odd mass was varied in the 
range 20 — 400 GeV. The upper bound is taken to ensure 
a proper separation between the light and heavy scales. 



as required by the EFT analysis. The very low mass 
range is expected to be severely constrained, but we defer 
the study of such region to future work. We turn next 
to a detailed description of the most important physical 
characteristics of the set of models in the scan, starting 
with the Higgs spectrum. 



III. MASSES OF LOW-ENERGY HIGGS 
BOSONS 

In this section we study the spectra of these models, 
analyzing the modifications with respect to the MSSM. 
Compared to the results already presented in 8], we in- 
clude the 1-loop supersymmetric corrections to the Higgs 
quartic couplings as given in [29] (a minor effect for the 
relatively low SUSY spectrum we have in mind), as well 
as the constraints coming from collider data (LEP and 
Tevatron) . 

In Fig. [1] we show the mass of the lightest CP-cvcn 
Higgs (h) as a function of m^, for both tan/3 = 2 
(left panel) and tan/3 = 20 (right panel). The green 
points represent models ruled out by LEP, while the 
magenta ones are excluded by current data from CDF 
and DO. We divide the remaining allowed models into 
two subsets. First, one has those models that will be 
probed at the Tevatron at 95 % C.L (red points), as- 
suming 10 fb~ per experiment and 50% efficiency im- 
provements [33] (see [37] for detailed projections in the 
MSSM context). These comprise the future reach of two 
search channels: h/H — >■ bb with the Higgs being pro- 
duced in association with electroweak gauge bosons, and 
gg — > h/H -^ W^W~ . Second, the blue points are those 
that will be out of the Tevatron reach under the previ- 
ous assumptions. For reference, we also show the MSSM 
(dashed) curve, assuming the same light SUSY spectrum. 
This color code will be employed in all the plots. 

The corrections to rrih due to the new physics are most 
important in the low tan/3 regime. Nevertheless, it is 
clear that they can also be relevant at large tan/3. The 
higher-dimension operators affect mh in such a way that 
it can easily be above the MSSM value. In the left plot, 
where tan /3 = 2, all the points lie above the MSSM curve; 
rrih can reach values as high as 250 GeV. 

Moreover, the left panel of Fig. [T] shows in a clear way 
how the Tevatron probes these models. For high enough 
values of ttia one distinguishes mostly uniformly single 
colored horizontal stripes. The magenta one, where nih ~ 
160 — 170 GeV, corresponds to h being excluded by the 
current Tevatron search in the WW channel '25]. Note 
that this range is slightly larger than the SM one (162 — 
166 GeV). This is due to the fact that, in our models, 
the gluon fusion cross section can be mildly enhanced 
with respect to the SM one. By the same token, one 
understands the presence of a few red points within the 
magenta stripe as those corresponding to models whose 
gluon fusion cross section is below the SM value. The 
two red stripes {rrih in the ranges 140 — 160 GeV and 
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FIG. 1. Lightest CP-even Higgs boson mass as a function of itia, for tan/3 = 2 (left panel) and tan/3 = 20 (right panel). We 
show the points excluded by LEP (green), excluded by current Tevatron data (magenta) and the region that will be probed by the 
Tevatron in the near future (red). The blue points are allowed by all the current experimental constraints. The dashed line is 
the MSSM result for the given SUSY spectrum. 



180 — 190 GeV) represent the future Tevatron reach of 
the h — >■ W^W~ channel. Notice also the presence of a 
thin stripe of red points, with vnh around 120 GeV, that 
extends along a wide range of rriA '. these models can be 
probed by the h -^ bb channel, that is effective only for 
relatively low values of m^. No points are excluded by 
the H ^f bb decay mode, since H is always much heavier 
than 120 GeV. 

The two blue stripes correspond to points where there 
is no reach from the Tevatron in the WW channel. This 
can be explained either by a low signal due to the reduced 
branching fraction into gauge bosons, or simply because 
the gg parton luminosity is not enough to produce such 
a heavy Higgs boson. Note however that in the high nrih 
blue region the ZZ — >■ Al channel becomes kinematically 
accessible, so that this Higgs could be observed in the gold 
plated four-lepton mode at the LHC. We will postpone 
further comments on this region to the next section. 

Regarding the LEP constraints, one sees that there are 
a few currently allowed (blue and red) points below the 
SM LEP bound of 114.4 GeV[3Sl|311- The nonexclusion 
is due to the fact that the coupling of h to the gauge 
bosons is reduced with respect to the SM value. However, 
all the red points below the LEP-bound can potentially 
be excluded in the H — )• WW channel. 

For the remaining points (m^ < 160 GeV, 
114.4 GeV < rrih < 170 GeV), the situation is more 
complex, and magenta, blue and red points coexist in 
this region. In particular, there is a region of allowed 
(blue) points with mh ~ 130 GeV — 140 GeV and rela- 
tively low mA- These points have suppressed branching 
fractions into both WW and bb, with AA being the dom- 
inant decay channel. 



In the case of tan/3 = 20, the deviations from the 
MSSM are far less dramatic. Ultimately, this is ex- 
plained by the fact that several higher-dimension oper- 
ators are tan/3 suppressed. However, m^ can reach val- 
ues as high as 160 GeV. In this case, since h is SM- 
like, the LEP bound is very strict, forcing rrih to be 
above ~ 110 GeV. Regarding the Tevatron searches, we 
see that there are two small and disjoint currently ex- 
cluded (magenta) regions. The region with rrih around 
160 GeV corresponds, as in the low tan/3 case, to exclu- 
sion based on the h — ^ WW decay mode. The second 
magenta region has lower values of ruh (114 — 130 GeV) 
and mA (100 — 135 GeV). This latter set of models are 
currently excluded by the inclusive tau search with 2.2 
fb~^, using the combination from CDF and DO [26 . This 
channel becomes important here, since the H/A -and in 
some cases the h- coupling to down- type fermions is tan /3 
enhanced. [? ] Turning to the red points (i.e. those within 
future Tevatron sensitivity), a closer inspection reveals 
that all of them can be excluded due to the decay modes 
of the lightest Higgs. In more detail, the ft, — >■ &6 channel 
probes points with mh below 126 GeV, while the rest are 
probed by the h — > WW search. This can be understood 
from the fact that, as in the MSSM, in the large tan/3 
limit H tends to be non SM-like. In contrast to the low 
tan/3 case, all the Tevatron allowed (blue) points corre- 
spond to somewhat heavy values oirriA (above 140 GeV). 

It is also interesting to study the relation between the 
CP-even Higgs masses. In Fig. [2] we show m,H as a func- 
tion of mil, for tan/3 = 2 (left panel) and tan/3 — 20 
(right panel) . For most of the points these masses are not 
correlated. For instance, if in the left plot one takes m^ 
in the 120 — 200 GeV range, then Mh can vary between 
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FIG. 2. tuh as a function of nih, for tan 13 — 2 (left panel) and tan/3 — 20 (right panel). We show the points excluded by 
LEP (green), current Tevatron data (magenta) and the region that will be proved by the Tevatron in the near future (red). The 
blue points are allowed by all the current experimental constraints. The dashed line is the MSSM result for the given SUSY 
spectrum. 



200 an(i 400 GeV. For tan/3 = 2, one has not only the 
(now vertical) stripes corresponding to exclusion due to 
h that we have found in Fig. [T] there are also horizontal 
stripes, corresponding to mn ranges where the Tevatron 
is excluding models by means of the H — )■ WW decay 
channel. This sheds some light into the region already 
mentioned in the discussion of Fig. [l] with m^ < 160 GeV 
and 114.4 GeV < ruh < 170 GeV. In this region both 
h and H can couple to WW, typically resulting in some 
suppression with respect to the SM for one or the other 
CP-even Higgs boson. This constitutes an interesting ex- 
ample of how the h and H signals can complement each 
other. The right panel confirms what we have antici- 
pated from our discussion of Fig. [T] in the large tan/3 
regime, uih tends to be heavy, and the decays of H are 
less restrictive than the ones from h, hence there are no 
horizontal stripes in this plot. As mentioned before, the 
h search channels give rise to all the red points. 

Finally, we show in Fig. [3] the masses of the heavy CP- 
even and charged Higgs bosons as a function of mj^. The 
deviations from the MSSM value are much less dramatic 
than for h. This is particularly true in the large rriA limit 
and for large tan/3. Nonetheless, in this region the con- 
tribution from the new physics effects to the masses is of 
0(10 GeV), which cannot be neglected. For low values of 
tan /3 (left plots) we see that in the region of blue points 
with low values of tua, both mH± and mn are above 
the MSSM value. Notice that this effect is more impor- 
tant for niH than for mjj± . As a direct consequence, in 
the low tan/3 regime, new exotic channels like H — ;> AA 
and H^ — > AW^ can be open, with large BRs, as we 
will see in the next section. This does not happen for 
tan/3 = 20 since, as stated before, there are no allowed 



(blue and red) points with rriA below 140 GeV and the 
mass splittings do not allow the previous decay modes. 

Having analyzed the modifications in the spectra due 
to the higher-dimension operators, we will devote the 
next section to study the collider phenomenology of these 
models. 



IV. BIVISSIVI COLLIDER PHENOMENOLOGY 

In this section we study the phenomenology of the 
BMSSM Higgs sector, including all of the effects and con- 
straints described in Section [Hi We consider the low and 
large tan/3 cases separately. 



A. Low tan/3 searches: general features 

We start with the low tan/3 regime, fixing tan/3 == 2. 
As we have described in the previous section, the main 
modification introduced by the higher-dimension opera- 
tors is to shift the Higgs spectrum with respect to the 
MSSM one. However, the couplings of the Higgs bosons 
also get corrected. The combination of these two effects 
can give rise to sizable modifications in both the Higgs 
production cross sections and the branching fractions. 

We compute the production cross sections in the fol- 
lowing way. For the Higgs-strahlung and vector boson 
fusion processes, we simply scale the corresponding SM 
cross section by the (square of) the Higgs- VK-VF coupling 
in our scenario, normalized to the SM coupling, i.e. by 
the effective coupling as defined in Eq. (IT]) (for all prac- 
tical purposes this ratio coincides with the normalized 
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FIG. 3. H (upper row) and H (lower row) masses as a function of tua, for tan/3 = 2 (left panels) and tan/3 = 20 (right 
panels). We show the points excluded by LEP (green), excluded by current Tevatron data (magenta) and the region that will be 
proved by the Tevatron in the near future (red). The blue points are allowed by all the current experimental constraints. The 
dashed line is the MSSM result for the given SUSY spectrum. 



Higgs-Z-Z coupling [? ]). For tlie gluon fusion cross 
section, we shall argue that the NLO K-factor in our sce- 
nario is expected to agree with the NLO K-factor in the 
SM within 20%. This implies that to this accuracy 



<y^lf[gg^h)^ T\^[h^gg) 



(2) 



since the ratio of cross sections equals the ratio of widths 
at leading order in as [SH HOI [JT] . The right-hand side 
of Eq. (pi) is computed using our modified version of 
HDECAY [27], which includes the tree- level expressions 
for masses and couplings in the presence of the higher- 
dimension operators. 



The K-factor in our scenario differs from the SM one in 
two respects. First, the contribution to the gluon fusion 
cross section from bottom loops cannot be neglected, spe- 
cially in the large tan/3 regime. Second, one has to con- 
sider the presence of a relatively light SUSY spectrum. 
We discuss separately these two effects. To assess the 
impact of the bottom loop we use the code HIGLU [H] , 
that includes both the LO and NLO results for both the 
SM and the MSSM [IT (but, at present, does not include 
SUSY particles in the loop), to compute the K- factors in 
these two models. We find that at low tan/3 and for 
a wide range of Higgs masses, the NLO K-factors for 
/i, H and A coincide within 5% with the SM NLO K- 



factor for a Higgs of the corresponding mass. At larger 
tan/3 (~ 30) the differences are larger, as expected, but 
still smaller than about 20%. We expect that the same 
will hold in our extended SUSY scenarios. The changes 
in the NLO K-factor due to relatively light sparticles in 
the loop, again in the MSSM context, were studied in 
[28) . where the effect was found to be less than 3% for 
tan/3 = 1.5. Therefore, we conclude that at low tan/3 
Eq. (pi) holds to an accuracy of better than 10%, and 
allows us to obtain a sufficiently precise estimate for the 
NLO gluon fusion cross section in our scenario. Note that 
this uncertainty is below the one obtained by comparing 
the NLO and NNLO/NNLL results in the SM calculation 
[43l - l48| . It is also important to note that the bulk of the 
effects of the light SUSY spectrum is taken into account 
in the LO cross section, and that these effects are fully 
implemented in HDECAY, which is used to compute the 
right-hand side of Eq. ^. This also includes radiative 
effects that correct the bottom Yukawa coupling, which 
can be important at large tan/^ [3T]. We discuss next a 
number of general features regarding the BMSSM Higgs 
signals. In both the SM and the MSSM, the dominant 
decay channel for a Higgs boson whose mass is greater 
than 140 GeV is into W pairs. Therefore, an important 
observable at a hadron collider is the production cross 
section times the branching fraction in the WW channel. 
In Fig. [4] we show this quantity for h and H, normalized 
to the SM result, as a function of the corresponding Higgs 
mass. In the left panel we clearly see the Tevatron ex- 
clusion in the h — >■ WW channel: the V-shaped magenta 
and red regions around mh '-^ 160 GeV correspond to 



the stripes that were already discussed in Section III (see 
Figs. Il]and[2]). We stress that the blue points with m^ 
above 180 GeV have a slightly enhanced WW signal com- 
pared to the SM. In turn, this mass range will be explored 
at the LHC via the h -^ ZZ — > 4/ channel (recall that, 
for all practical purposes, the CP-even Higgs normalized 
couplings to WW and ZZ are the same, hence the plot 
can be directly applied to the ZZ channel). Thus, an 
enhanced signal in this region is an interesting feature: 
for these points, the Higgs cannot escape detection. For 
the blue points with WW signal reduced by a factor of 
10 or more {mh < 160 GeV), one may have to rely on 
other search channels. 

Note that this figure exhibits currently allowed (blue 
and red) points with rrih below the LEP bound of 
114.4 GeV. These correspond to models where the cou- 
pling to gauge bosons is below the SM value. We also 
notice a group of red points whose signal is around the 
SM value, and with a mass slightly above the LEP bound 
(114.4 GeV < rrih < 120 GeV): these are within the 
Tevatron reach in the h ^>- bb channel, assuming an accu- 
mulated luminosity of 10 fb^ per experiment and a 50% 
efficiency improvement in this channel. 

Turning our attention to the right panel of Fig. [41 we 
see that in the case of H it is hard to differentiate regions 
where a single color is predominant, as was possible in 
the left panel. We can identify a mostly green region 



with rriH above the LEP bound (and above the MSSM 
curve). These points are excluded by the LEP bound 
on m^ rather than on mn , and serve as a reminder that 
the constraints may come from observables not related to 
those shown in a given plot. This is not to say that there 
are no points where the exclusion is through H directly 
instead of h: for instance, the magenta and red points in 
the upper left side of the plot correspond to the V-shape 
exclusion from H — > WW at the Tevatron. This is the 
only region where the signal is enhanced with respect to 
both the SM and the MSSM. These correspond to models 
where H is SM-like, while h decays mainly into bb and 

TT. 

Aside from the WW channel, there are other impor- 
tant decay modes for light Higgs bosons, in particular 
bb, TT and 77. In the first case, the huge QCD back- 
grounds render this channel very difficult to measure at a 
hadron collider. This does not mean, however, that this 
decay mode is completely useless. For instance, in the 
Higgstrahlung process, qq -^ Z* jW* — ^ ZjW + Higgs, 
the gauge boson can be fully reconstructed from its de- 
cay modes, and then Higgs — >■ 65 becomes a feasible op- 
tion. Another example to search for a SM-like Higgs 
boson decaying into 6& is the Higgs associated produc- 
tion together with a top quark pair. This has the prob- 
lem of being quite challenging at hadron colliders. The 
di-photon channel, on the other hand, constitutes the 
most promising decay channel for a relatively light SM- 
like Higgs at the LHC since, in spite of its tiny BR of 
0(10^"^), an excellent energy resolution can be achieved 
and the background is under good experimental control. 
The other important search channel at the LHC in the 
low Higgs mass range is the vector boson fusion with the 
subsequent decay of the Higgs into a tt pair. 

In Fig. p^ we show the branching fraction of h into hb 
(left plot) and 77 (right plot), for tan/3 = 2. Notice 
that the hb channel can be suppressed with respect to 
the SM one, as in the blue points with masses in the 
120 - 150 GeV range and BR(/i -^ bb) < 10"^ This is 
an interesting feature, since it can lead to enhancements 
in other search channels. One can also see currently al- 
lowed (blue and red) points with BRs into bb above the 
SM curve: those have a reduced BR into W's, as we have 
previously identified in Fig. [4] In the case of H (not 
shown here), the BR into bb is typically higher than the 
SM value. With respect to the MSSM, we find that there 
is no definite tendency: BR(_ff — > bb) can be either in- 
creased or suppressed by an order of magnitude. It is 
worth mentioning that the branching fraction in the tt 
channel follows closely the bb behavior. This is as ex- 
pected, since the extended Higgs sectors under consider- 
ation do not distinguish between the down-type fermions, 
in the sense that the Yukawa coupling normalized to the 
SM value is the same for bottoms and taus, while differ- 
ences due to the SUSY QCD and top Yukawa interac- 
tions, that arise at loop level, are not significant at small 
tan/3. 

Turning our attention to the right panel of Fig. [5l 
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FIG. 4. Production cross section by gluon fusion times branching ratio into W boson pairs, for h (left panel) and H (right 
panel), normalized to the SM result, as a function of the corresponding Higgs mass, for tan/3 = 2. We show the points excluded 
by LEP (green), excluded by current Tevatron data (magenta) and the region that will be probed by the Tevatron in the near 
future (red). The blue points are allowed by all the current experimental constraints. The dashed line corresponds to the MSSM 
prediction for the given SUSY spectrum. 
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FIG. 5. Branching fractions for h ^f bb (left panel) and ft — > 77 (right panel) for tan/3 = 2. We show the points excluded by 
LEP (green), excluded by current Tevatron data (magenta) and the region that will be probed by the Tevatron in the near future 
(red). The blue points are allowed by all the current experimental constraints. The solid (dashed) line corresponds to the SM 
(MSSM) result. 



we see that most models present a suppressed branching 
fraction in the di-photon channel. However, it is worth 
noticing the group of points above the SM curve, where 
an enhancement of up to a factor of 2 can be achieved. 

Decays of the CP-even Higgs bosons into pairs of A 
bosons can become the dominant decay mode. Such a 
scenario has been previously considered in the literature 
(see, for instance, [IS] for a model-independent analysis, 
and Uni for NMSSM studies). In Fig. [6] we show the 



branching fraction of h and H into AA, for tan/3 = 2. 
The left panel shows that the branching fraction in this 
channel can reach 0(1) values, thus becoming the most 
relevant decay mode of h. The Tevatron allowed (blue) 
points in this figure present a reduced branching fraction 
in both the bb and the WW channel, and were already 
mentioned in the context of Fig. [5] In the case of H 
(right panel), the branching fractions vary considerably, 
but the AA channel may still become the primary decay 
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FIG. 6. Branching fractions for h — >■ AA (left plot) and H -> AA (right plot), for tan/3 = 2. We show the points excluded by 
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for the given SUSY spectrum. 



mode in some models. 

Tlie main modification to tlie decay plienomenology of 
A and H^ witli respect to the MSSM is due to tlie shift in 
the overall Higgs spectrum. The channels that change the 
most are those that involve a Higgs decaying into either 
a pair of Higgs bosons, or a Higgs boson plus a gauge 
boson. As an example of the latter, we take the decay of 
H^ into a W boson and a neutral Higgs. In the MSSM, 
since h tends to be rather light, one has that H"^ — >• hW^ 
can be an important decay channel. On the other hand, 
since niA and m.H± tend to be rather degenerate, one 
finds that H^ — > AW^ is generally highly suppressed. 
On the contrary, in the context of the BMSSM, one can 
find points where the mass hierarchy suffers an inversion, 
i.e. rriji can be well below tti/j and split from mH± . In this 
case, one finds that A and h interchange their roles with 
respect to the above described situation in the MSSM, as 
can be seen in Fig. [7] The left panel shows points where 
the AW^ channel has a BR greater than 0.1, while the 
right panel shows that the hW^ decay mode is highly 
suppressed. In this case, the process H"^ -^ AW^ — 7> 
bbW^ can give rise to an interesting signal at the LHC, 
possibly allowing the discovery of two nonstandard Higgs 
bosons. Note that we have already encountered another 
example of an inversion between A and h in the context 
of Fig. |6] the potentially open MSSM channels A ^ hh 
and H ^ hh are replaced by h -^ AA and H — > AA in 
the BMSSM context. 

Finally, turning our attention to A, we have found that 
the A — >■ hZ decay channel is significantly reduced with 
respect to the MSSM value ^ 0.3 for values of m^ below 
250 GeV, due to the shift in nih that disfavors this decay 



mode. [? ] This reduction brings an enhancement in both 
the bb and tt channels. As in the MSSM, the former is 
the dominant decay channel below the ti threshold, and 
the latter stays almost constant at about 10%. 



B. Low^ tari/3 searches: benchmark points 

Up to this point we have analyzed each observable 
almost independently of the others. We would like to 
understand, however, how the different features that we 
have singled out are correlated with each other. We shall 
consider benchmark scenarios currently allowed by LEP 
and Tevatron data and explore two possibilities: a) mod- 
els that can be probed at 95 % C.L at the Tevatron in the 
near future, from now on referred to as Tevatron covered 
(red) points, and b) models that are beyond the expected 
Tevatron reach and will be explored at the LHC, from 
now on referred to as Tevatron uncovered (blue) points. 
We will also indicate the importance of the 1/M^ effects, 
as measured by their impact on nih- 



1. Scenarios within the Tevatron reach 

The Tevatron covered models can be divided into three 
subsets, according to which channel can exclude the 
point: h -^ bb, h ^ WW and H -^ WW. It is in- 
teresting to ask whether a given model can be probed 
by more than one channel at the Tevatron. We find, 
however, that the previous subsets are disjoint. The dis- 
jointness between the subsets probed by /i — > 66 and 
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hW (right panel), for tan ,(3 = 2. The dashed line corresponds to the MSSM 



h — >■ WW can be understood in terms of the relevant 
mass ranges, since the hb search is most sensitive to the 
rrih ^120 GeV range, while the di-boson channel probes 
the region 165 ± 20 GeV. In principle, they do not have 
to be mutually exclusive, but one would need an en- 
hancement of 3.4 over the SM in the h — > WW signal 
in order to probe a 120 GeV Higgs in this channel, [? 
] which is not achievable within these models: the in- 
crease in the gg —¥ h -^ WW cross section is always 
below a factor of 2. In the case of h and H decaying 
into W bosons, even if both of them are in the favorable 
mass region (-- 150 - 170 GeV), the MSSM sum rule 
9hww + 9hww = 1 is valid within 5 % accuracy, and 
it is not possible for h and H to have large enough cou- 
plings to W's for both signals to simultaneously be within 
the Tevatron reach. The subsets probed by the ft, — >■ 66 
and H — > WW searches are disjoint because both pro- 
cesses require a sizable coupling of the Higgs to W^'s (for 
production and decay, respectively), but this does not 
happen when ttt,^ and niff are sufficiently different, as 
would be required for simultaneous searches in these two 
channels. 

Point A: MSSM-like scenarios 

We start our analysis with the points that can be 
probed via the h — > bb decay mode. We find that these 
do not differ greatly from the decoupling limit of the 
MSSM with rather heavy sparticles (~ a few TeV). In 
this case, the observation of a light SUSY spectrum (in 
the few hundred GeV range) would be the smoking gun 
of BMSSM physics, since such a light SUSY spectrum 
would be in conflict with the LEP limits on the MSSM 
lightest CP-even Higgs boson. We illustrate the main fea- 
tures of this subset by showing point A [? ] in Table llj 
where we include the mass spectrum and the branching 



POINT A 



rriA (GoV) 


mfe (GoV) 


mn (GcV) 


m^± (GcV) 


239 


118 


246 


245 


2 

9hww 


2 

Shww 


9hgg 


ffffgg 


0.992 


0.008 


1.06 


0.55 


channel 


BMSSM (SM) 


channel 


BMSSM (SM) 


h^ bb 


0.78 (0.73) 


h -i. WW 


0.08 (0.11) 


h^ TT 


0.08 (0.08) 


h -^ 77/10"''' 


1.42 (2.30) 


H ^ bb 


0.15 


H -^ WW 


0.22 


H ^ ZZ 


0.11 


H -^ hh 


0.50 


A^ bb 


0.89 


H+ -^ tb 


0.99 


A^Tt 


0.08 


A^ Zh 


0.24 



TABLE I. Masses and branching fractions in the BMSSM 
(and m the SM for h) for point A. We only show the main 



decay modes. 
Eq. 0). 



The effective couplings g^x were defined in 



fractions of the most important decay channels for each 
Higgs boson. For reference, in the case of h we also in- 
dicate between parentheses the SM values. We also note 
that for this point, the 1/M^ operators contribute about 
25 GeV to ran. 

Generically, the branching fractions of h do not devi- 
ate much from the SM ones. One finds a small increase 
in the h ^^ bb channel, while BR(ft — > 77) is slightly sup- 
pressed with respect to the SM (by at most a factor of 
3). Since gf^^/y^ and g^ are close to one, the produc- 
tion cross sections by Higgs-strahlung and gluon fusion 
are, for all practical purposes, the same as in the SM. 
Thus, the change in the signal is given by the ratio of 
the branching fractions in our scenario to those in the 
SM. For point A, the production rate in ft — >■ 66 is 6% 
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POINT B 



POINT C 



■niA (GoV) 


mn (GeV) 


mn (GeV) 


ni„± (GeV) 


101 


129 


141 


135 


9h\v\v 


2 

9hww 


2 
9hqa 


Bhciq 


0.8 


0.2 


1.72 


0.06 


channel 


BMSSM (SM) 


ehannel 


BMSSM (SM) 


h^bb 

h -^ WW 

h — >■ jets 

H -^ bb/rf 

A -^ bb/rf 


0.01 (0.56) 
0.63 (0.28) 
0.26 (0.06) 
0.84 / 0.09 
0.89 / 0.09 


h^ Tf 

h^ ZZ 

/i^77/10-^ 

H -> WW 


0.001 (0.06) 
0.08 (0.04) 
3.97 (2.38) 

0.05 

0.87 



TABLE II. Masses and branching fractions in the BMSSM 
(and in the SM for h) for point B. 



above the SM result. In this case, the Tevatron could 
claim a hint on a SM Higgs boson, while at the LHC 
the direct detection of h would proceed in the di-photon 
channel, since the gg ^ h ^ ^^ cross section is 0.65 of 
the SM value. Some of the remaining Higgs bosons may 
also be observed. For H and A, the H ^ hh ^ jjbb and 
A — > Zh — >• llbb searches provide the best prospects for 
discovery [52]. For a charged Higgs with a mass above 
rrit, the ATLAS update of 2009 53] found that the tb 
channel is rather challenging and that the low tan/3 re- 
gion cannot be covered. 

Point B: Light Higgs spectra 

We turn now our attention to the models that can 
be excluded at the Tevatron by the h — >■ WW channel. 
Those can be further split into two categories, according 
to whether ruh is high (> 170 GeV) or low (< 160 GeV), 
corresponding to the two red stripes defined in the con- 
text of Fig. [T] As a general feature of the lower red 
stripe, the branching fraction of h into bb can be sizably 
reduced with respect to the SM, as we pointed out in 
the left panel of Fig. [5] This implies that the remain- 
ing channels are enhanced, which is interesting for the 
h — )■ WW and /i — J- 77 decay modes. We present as an 
example point B [? ] in Table Ih] Here, one sees that 
the Higgs spectrum is relatively light. H is the heaviest 
Higgs, while h is lighter than H^, but heavier than A. It 
turns out that for this point, the 1/M^ effects result in a 
slight net reduction of rrih by a couple of GeV. 

Since h is SM-like, we give in parentheses the corre- 
sponding branching fractions in the SM. Here we clearly 
observe that h presents an increase in the gluon fusion 
cross section, and in the branching fractions into pho- 
tons and W bosons, accompanied by a sizable reduction 
in the down-type fermion decay modes. Note that the 
53 — >■ /i — >■ 77 and gg —>■ h —i' WW signals are larger 
than in the SM by factors of 2.86 and 3.82 respectively, 
which would facilitate the search of h at the LHC as well. 
H decays mainly into bottoms and taus, and its produc- 
tion cross section by gluon fusion is strongly reduced with 
respect to the SM case. The most promising discovery 
channel at the LHC would be qqH — > qqrf, where the 



■niA (GeV) 


nih (GeV) 


niH (GeV) 


ni„± (GoV) 


135 


174 


186 


164 


9h\v\v 


9h\v\v 


2 

9hqq 


9Hqq 


0.11 


0.89 


1.05 


0.65 


channel 


BMSSM (SM) 


channel 


BMSSM (SM) 


h^ bb 

H -^ WW 

A^ bb 


0.12 (0.01) 
0.81 (0.82) 

0.90 

0.59 


ft-i- WW 
H ^ ZZ 

A^ Tf 

H+ -J- tb 


0.84 (0.96) 
0.17 (0.17) 

0.10 

0.38 



TABLE III. Masses and branching fractions in the BMSSM 
(and in the SM for h and H) for point C. 



signal is reduced with respect to the SM by a factor of 
two. The CP-odd A decays as in the MSSM, while for 
the charged Higgs the tv^ channel is the dominant one. 

Point C: The heavy CP-even H as the SM-like Higgs 

In the high mh region that can be probed at the Teva- 
tron in the h — )■ WW channel one finds an unusual SUSY 
spectrum. Typically, one runs into the previously men- 
tioned inversions between h and A. Moreover, h can also 
be heavier than the charged Higgs, which is a feature that 
is not present in the region where m^ is below 150 GeV. 
We illustrate this with point C [? ] in Table HI For this 
point, the 1/il/^ operators contribute about 30 GeV to 
ruh- 



Here we see that the two CP-even Higgs bosons have 
masses well above the maximum value for rUh that can 
be obtained in the nifi max scenario within the MSSM 
context. In this case, it makes sense to compare both h 
and H with the SM Higgs. For this particular point, h 
has not been yet excluded by the Tevatron search since 
it is not SM-like and its branching fraction into WW is 
somewhat suppressed. H can be discovered by the LHC 
in the ZZ — >■ 4/ mode, since the signal normalized to 
the SM value, is 0.65. We note that here both rrih and 
mf{ are near the region where the WW channel opens 
up leading to a suppression in the sensitivity of the ZZ 
search mode at the LHC. We recall that such a heavy 
SM-like H is not a feature of the MSSM, being a unique 
characteristic of the BMSSM Higgs sector. The CP-odd 
A decays almost entirely to bottom and tau pairs, while 
the charged Higgs has sizable decays into both the ri'r 
and tb channels. 

The last subset of the Tevatron covered points corre- 
sponds to those than can be probed by the H — >■ WW 
search, for which the gg —>■ H —i' WW signal goes be- 
tween 0.4 — 4 times the SM value. In such scenarios, h 
and A decay mostly into bottoms and taus. In some cases 
the h ^ "f"f OT qqH — ^ qqrr signal might be observable 
at the LHC. The charged Higgs is relatively light (al- 
ways below 200 GeV) and will decay almost 100% of the 
time into ri'r for masses below 160 GeV, and in tb for 
the remaining points. The H^ —>■ hW^ decay mode is 
closed for kinematical reasons, as we already know from 
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POINT D 



POINT E 



■niA (GcV) 


mu (GcV) 


rriH (GcV) 


m„± (GcV) 


184 


204 


234 


203 


2 
9hWW 


2 

9/f WW 


9hqq 


2 
9Hqq 


0.3 


0.7 


1.39 


0.36 


channel 


BMSSM (SM) 


channel 


BMSSM (SM) 


/i-)- WW 
ff -i- WW 

A^bb 


0.73 (0.72) 

0.70 (0.71) 

0.87 


h^ ZZ 
H ^ ZZ 

H+ -i. tb 


0.25 (0.27) 

0.29 (0.29) 

0.99 



TABLE IV. Masses and branching fractions in the BMSSM 
(and in the SM for h and H), for point D. 



Fig. [7| In addition, when m,A is light, the h -> A A and 
H^ -^ AW^ channels might become important. Since 
this also happens with the Tevatron uncovered points, we 
will defer further comments and the study of a suitable 
benchmark point for the next subsection. 



2. LHC searches 

Regarding the Tevatron uncovered points, we can also 
split them into two disjoint subsets, corresponding to 
each of the blue stripes in Figs. [T] or [2} we will refer to 
them as low mass {ruh < 140 GeV) and high mass (above 
190 GeV) regions. In the high ruh case, one can make 
a further distinction according to whether tha is below 
or above 160 GeV. Again, we illustrate the possibilities 
with a few benchmark points. 

Point D: Two peaks in the ZZ — > 4/ signal 

General features of the high ?n/i, high jtia case are: 
increased gg —)' h cross section with respect to the SM, 
and negligible (below 2%) changes in the h — >■ WW/ZZ 
decay modes. Regarding H, one has that the signal in the 
gg —>■ H —i' WW/ZZ channel is always suppressed with 
respect t o th e SM. As an example, we show point D [? ] 
Here, the 1/Af^ operators contribute about 



IV 



in Table 

40 GeV to rrih- Given the features of this point, it again 
makes sense to compare both CP-even Higgs bosons with 
the SM. 

The rise in the mass of h automatically closes the 
H — > hh, A — > hZ and H^ — >■ hW^ decay modes, which 
could be important in the MSSM case. This picture 
can suffer some alterations if h stays around 200 GeV, 
while the rest of the Higgs bosons attain values around 
400 GeV, since this will open not only the previously 
mentioned channels, but possibly also decays into sparti- 
cles. In such a case, one would run into a sort of MSSM 
decoupling limit, but with a mass for the lightest Higgs 
which is unattainable within the MSSM. Concentrating 
on point D, we emphasize that both h and H couple 
in a sizable way to the electroweak gauge bosons, and 
thus the measurement of both couplings will permit a de- 
tailed study of the EWSB mechanism, as it arises from a 
2HDM. One possibility to discover these CP-even Higgs 



niA (GeV) 


rah (GeV) 


ran (GeV) 


m„± (GcV) 


134 


181 


205 


165 


2 

9h \v \v 


2 

9hv/v/ 


9hqq 


2 
9Hqq 


0.03 


0.95 


0.79 


0.99 


channel 


BMSSM (SM) 


channel 


BMSSM (SM) 


h^bb 
/i-s- WW 
H -^ WW 

A^bb 
H+ -^ tb 


0.23 (0.005) 
0.68 (0.92) 
0.72 (0.73) 

0.89 

0.57 


/i — >■ rr 
h^ ZZ 
H ^ ZZ 

A^ Tf 
H+ -^ TV^ 


0.03 (0.0005) 
0.04 (0.07) 
0.27 (0.27) 

0.10 

0.40 



TABLE V. Masses and branching fractions (and in the SM 
for h and H) for point E. 



bosons would be to search for two isolated peaks in the 
ZZ -^ 4/ golden mode. Notice that the branching frac- 
tions of both h and H are very close to their SM coun- 
terparts, while there is a difference in the gluon fusion 
production cross section. Since the gg ^ H ^ ZZ sig- 
nal can be sizably suppressed with respect to the SM, 
the direct detection of H in this channel might not be 
feasible. Nevertheless, a large number of models similar 
to point D would in fact present two clear peaks in the 
di-lepton invariant mass distribution. We also note that 
these points have a gg —i' h —)' WW/ZZ signal that can 
be 20-70 % larger than the SM value. 

Point E: Non-SM-like Higgs with a clear di-hoson signal 
One interesting example of a point where m^ is still 
high, but ruA is below 160 GeV is given in point E [? ] , 
shown in Table |VJ We see that this point has an unusual 
Higgs hierarchy, since here h is heavier than both A and 
H^ (in this respect similar to point C). The signal for 
h -^ WW, normalized to the SM, is 0.54. The CP-even 
H will be discovered first and will appear to be the SM 
Higgs, since the signal is very close to the SM one. Soon 
after, h will be found in both the ZZ and WW channels, 
thus providing a clear evidence of new physics. Notice 
that here the coupling of h to gauge bosons is extremely 
small, but due to kinematics it still decays preferentially 
into gauge bosons. We stress again that this is a unique 
characteristic of the BMSSM Higgs sector in the low tan f3 
regime, since in the MSSM a similar behavior can only 
occur for H. For this point, the 1/M^ operators con- 
tribute about 30 GeV to to/j. 

Point F: Multi-Higgs decay chains 

For the points outside the Tevatron reach with tti^ be- 
low 140 GeV, the most remarkable feature is the pos- 
sibility of having the channels h -^ AA and H -^ AA 
kinematically open. For the points where these chan- 
nels are closed, the situation is not as interesting, so 
we will focus on the first scenario. As an example, we 
show point F ['i 



in Table VI 



where these channels 
are the dominant decay modes of both h and H. For 
this point, the 1/M^ operators contribute about 10 GeV 
to ruh- Focusing on the AA channel, the possible final 
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POINT F 



niA (GoV) 


mn (GcV) 


niH (GoV) 


ni„± (GeV) 


64 


135 


155 


125 


2 

9hww 


2 
5/f H^H' 


3hgq 


2 


0.002 


0.991 


0.65 


1.17 


channel 


BMSSM 


channel 


BMSSM 


h^bb 

H -^ WW 

H ^ bb 


0.15 
0.12 
0.02 
0.56 


h^ AA 

H ^ AA 

A^bb 

H"^ -j. w^ + A 


0.84 
0.84 
0.92 
0.40 



TABLE VI. Masses and branching fractions in the BMSSM 
for point F. 



states for h and H are 6666, 66rf and tttt jlQHoT] . The 
first one is very challenging due to the enormous QCD 
background, while the third one suffers from a reduced 
signal [BR{A — >• rf) ~ 10%]. This leaves the 66rf chan- 
nel as the most promising one. For the case of H, one 
may also look at the gg ^ H ^ WW channel, whose sig- 
nal is 1/5 of the SM value, and could be discovered with 
about 100 fb^^ [52] • For the charged Higgs, the dominant 



decay mode is rvr- Notice also that H^ 



AW^ 



have a sizable branching fraction, offering the possibility 
to discover both A and H^ simultaneously in this decay 
mode. 



C. Large tan/3 searches: general features 

In this subsection we present our analysis for the large 
tan/3 regime, fixing tan/3 = 20. As shown in Section [iTTl 
the changes in the spectrum with respect to the MSSM 
are less important than in the low tan/3 case. 

We use Eq. ^ to estimate the gluon fusion production 
cross section at NLO in a^ . Although the impact of the 
bottom loop in the A'-factor is more important for larger 
tan/3, the NLO K-factor in our model is still expected to 
be within 20% of the NLO SM K-factor, as discussed at 
the beginning of Section p^V A[ Furthermore, as shown in 
[55], the effects on the K factor due to a light sparticlc 
spectrum like the one we are considering are negligible at 
large tan/3. Hence, we conclude that simply computing 
the right-hand side of Eq. ^ allows us to obtain the 
NLO gluon fusion production cross section within 20% 
accuracy even at large tan/3. In this regime, production 
in association with a 66 pair can become important, and 
can be obtained in our model from existing results by a 
simple rescaling with the effective coupling (?5f,(,, where 
<j>^h,H,A [seeEq. Q]. 

We show in Fig. [8| the effective coupling of h to down- 
type fermions, .9?ie, for both large (left panel) and small 
(right panel) tan/3. At large tan/3, one sees that the 
currently allowed models (blue and red points) have a 
gluon fusion production cross section which ranges from 
0.7 — 5 times the SM value. The most striking feature 
is that the coupling to bottom pairs can be strongly 



suppressed. For large tan/3 this happens for a value 
of (T{gg — )■ h)/SM of around 1.4. These models have 
110 GeV < mil < 150 GeV, where the decays into 66 of 
the SM Higgs are important. The suppression in g^.j can 
also be observed at low tan/3 (right panel). In this case, 
the associated values of rrih are in the somewhat higher 
range from 120 GeV to 250 GeV, with the strongest sup- 
pressions occurring for ruh > 150 GeV. 

The suppression in the coupling to down-type fermions 
is somewhat reminiscent of the small aoff scenario [54| 
I55| . but there are important differences. In the small 
aeff scenario the 5^^^ coupling is suppressed as a result 
of a cancellation between the tree-level and one-loop con- 
tributions. This can happen at large tan/3, where the ra- 
diative effect is enhanced at the same time that the tree- 
level contribution is somewhat suppressed, thus allowing 
for a cancellation. Besides large tan/3, sizable values of 
jiAt/MgjjgY S'i'G necessary, and the cancellation is found 
to happen only for certain values of tua (below or of or- 
der 200 GeV) that are highly correlated with tan/3 |54j . 
In contrast, the suppression we find occurs as a result 
of a cancellation between the tree-level MSSM contri- 
bution and those due to the higher-dimension operators 
(we have checked that the picture remains unchanged 
by turning off all loop effects). Most importantly, the 
fact that the suppression occurs at tree-level implies that 
the couplings to bottom and tau pairs are simultaneously 
(and strongly) suppressed. This does not tend to happen 
in the small Ucs scenario, since the radiative enhance- 
ments for bottoms and taus happen in different regions 
of parameter space. Also, in spite of the large number 
of parameters, there is a clear correlation between the 
g^,T suppression and the <7{gg — >■ h) enhancement. The 
increase in the gluon fusion cross section is due to the 
destructive interference of the bottom loop in the gluon 
fusion cross section, and also to light SUSY particles run- 
ning in the loop. 

The suppression of the down-type fermion channels im- 
plies a general enhancement in the branching fractions of 
the remaining channels. The most interesting enhance- 
ments are those in the gauge boson channels: WW, ZZ 
and 77. In the left panel of Fig. [9|we show the branching 
fraction into WW at large tan/3, again as a function of 
the gluon fusion production cross section normalized to 
the SM value, which can be compared to the left panel 
of Fig. [sj We see that the region where the hbb cou- 
pling is suppressed is exactly where the WW branching 
fraction is greatly enhanced, and leads to an interesting 
Tevatron sensitivity in the W channel over a wide range 
of rrih ■ The left panel of Fig. [9| clearly exhibits how the 
Tevatron covered (red) points arise. The upper red re- 
gion corresponds to those models within Tevatron reach 
in the h —>■ WW search, while the two lower red regions 
contain only points that can be probed in the /i — >■ 66 
channel. These latter models always have a branching 
fraction into WW below about 20%. 

In the right panel of Fig. [9| we show the 517 — >■ /i — )• 77 
cross section, normalized to the SM. We see that the di- 
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tan /S = 20, M = 1 TeV, ii = m, = 200 GeV, Msusy = 300 GcV, A, =Ai,=0 
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FIG. 8. Effective coupling g^^^^ for tan/3 = 20 (left panel), and for tan/3 — 2 (right panel), as a function of the gluon fusion 
production cross section over the SM value. We show the points excluded by LEP (green), excluded by current Tevatron data 
(magenta) and the region that will be probed by the Tevatron in the near future (red). The blue points are allowed by all the 
current and near future experimental constraints. The dashed line corresponds to the MSSM result for the given SUSY spectrum. 
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(green), excluded by current Tevatron data (magenta) and the region that will be proved by the Tevatron in the near future (red). 
The blue points are allowed by all the current experimental constraints. The dashed line corresponds to the MSSM result for 
the given SUSY spectrum. 



photon signal can be increased with respect to the SM 
one by up to a factor of 10. This strong enhancement is a 
direct resuh of the decreased branching fraction into 66, 
together with the enhancement in the gluon production 
cross section discussed above. The points with enhanced 
signal in the diphoton channel correspond to values of rrih 
between 110 GcV and 130 GeV. It is interesting to com- 
pare to the latest available diphoton analysis from CDF 



(Sg and DO [57]. The CDF analysis, performed with 5.4 
fb~^ of data, quotes an observed hmit of 18.7-25.9 for 
the di-photon cross section normalized to the SM. The 
DO analysis, with 4.2 fb^^, gives a corresponding limit 
of 11.9-28.3 . [? ] As a result, the enhancement in the 
di-photon signal we find can be interesting at the Teva- 
tron, and of course it would be spectacular at the LHC. 
One should also notice that for models with enhanced 
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BR{h -^ bb), the signal into photons can be reduced by 
up to a factor of 10. 

The gg -^ h ^f WW signal (not shown here) presents 
the same behavior as the 77 one. This can be easily 
understood as follows. In the SM, the ft, — >■ 77 decay 
mode proceeds via W and top loops, the former giving 
the dominant effect. In our currently allowed (blue and 
red) points, the coupling of h to tops and W's is very 
close to the SM value (the differences are below 2%). 
Although gy^f^i can be enhanced by a factor of 10, the 
bottom loop is still a small contribution to the /i — >■ 77 
process. Therefore, the partial widths T(h — )■ WW) and 
T{h — )■ 77) in our model are very close to the SM ones, 
and the changes in the branching ratios of each channel 
are common and strictly due to the variation of BR{h — >■ 
bb) with respect to the SM. Therefore, enhancements in 
the WW/ZZ channels can also be interestingly large. 

With respect to the remaining Higgs bosons, the situa- 
tion resembles the large tan/3 regime of the MSSM. Both 
H and A decay mainly into bottoms and taus, while the 
charged Higgs goes to either r^'^- or tb depending on its 
mass. It is also possible for a heavy Higgs to decay into 
the lightest one: BR(iy -^ hh) can reach 30%, while both 
BR{A -^ hZ) and BR(iy± -^ hW^) can reach 10%, pro- 
vided the decaying Higgs boson mass is above 200 GeV. 
For this mass range, the decay mode into sparticles can 
also be important, if kinematically allowed. 



D. Large tan/? searches: benchmark points 

Having described the main differences of the large tan- 
gent beta regime with respect to the MSSM, we show a 
selected sample of benchmark points. 



1. Scenarios within Tevatron reach 

We start with the points covered in the near future by 
the Tevatron via the h ^ bb search. From Fig. [o] one 
sees that for these points (lowest red region in the left 
panel), the signal into 77 (and WW) can be enhanced 
by at most a factor of 2. Since such enhancement fac- 
tors can also be obtained within the MSSM (for instance 
with sparticle masses around 500 GeV) , we will not show 
a benchmark point here, but will briefly comment on the 
main characteristics of these type of models. The Teva- 
tron could claim a hint in the bb channel, while at the 
LHC the signals into 77 and rf are enhanced by up to 
a factor of 2 with respect to the SM, thus allowing for 
a discovery using these decays modes. Regarding the 
remaining Higgs bosons, one sees that both H and A 
decay mainly into bottoms and taus. The giuon fusion 
production cross section for H and A is around 80% of 
the SM-value, while the bbh production becomes an im- 
portant mechanism due to the large tan/3 enhancement 
(the CP-even H has highly suppressed couplings to W's 
and Z's). Thus a discovery in the H/Att search may be 
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BMSSM (SM) 


h~¥ bb 
h^ ZZ 
H ^ bb 
A^ bb 


0.43 (0.20) 
0.08 (0.05) 

0.75 

0.84 

0.21 


h -s- WW 

H ^ Tf 

A^ Tf 

ff± -^ tb 


0.07 (0.02) 
0.41 (0.66) 

0.13 

0.14 

0.75 



TABLE VII. Masses and branching fractions in the BMSSM 
(and in the SM for h) for point G. 



feasible [52]. We note also that A and H can be very 
close in mass, so that the two states cannot be disentan- 
gled at the LHC, but rather the signals have to be added 
up. Because of a sizable branching fraction into the rur 
channel, the charged Higgs can be within LHC reach, 
even for m//± > nit [IS]- The way to distinguish such a 
situation from the MSSM will be through the observation 
of relatively light superparticles. 

As we mentioned before, an interesting possibility is 
to have sizable branching fractions for the decay modes 
H -^ hh, A ^ hZ and H"^ -J> hW"^. This requires 
heavy Higgs bosons with a mass above 250 GeV. De- 
pending on the details of the SUSY spectrum, also de- 
cays into sparticles may be open. We have found that the 
branching fractions in these multi-Higgs channels are be- 
low 10% in most cases, and that one would still have both 
A and H decaying sizably into down-type fermions, with 
H^ decaying preferably into tb but with a non-negligible 
branching fraction into tVt due to the large tan/3 en- 
hancement. Provided that niH > 300 GeV, the H — ^ hh 
branching fraction can reach values of up to 20 — 30%, 
which is interesting since it allows for the potential ob- 
servation of several Higgs states. 

Point G: SM-like Higgs heavier than the MSSM upper 
bound 

We turn now to the models covered at the Tevatron 
via the h — )■ W W sea rch (upper red region in Figs. |9|. 

an example where ruh is above 



We show in Table VII 



the maximum attainable value in the m^ max scenario of 
the MSSM, with sparticles at the TeV scale. The 1/AP 
operators contribute about 45 GeV to m/j. We notice 
that in this point the enhancement in the gluon fusion 
production cross section [1.64 x <TsM{gg -^ h)] is com- 
pensated by the reduction in the WW branching fraction 
with respect to the SM (0.41/0.66), thus resulting in a 
gg ^ h ^ WW signal close to the SM one. As can 
be seen in the right panel of Fig. [o] (interpreted for the 
WW channel) this is a general feature of the Tevatron 
covered points in the higher range of mh ■ The nonstan- 
dard neutral Higgs bosons, H and A, can be detected in 
the TT (or /x/i) channels, as is well known for the large 
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POINT H 
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h^ bb 
h — f jots 
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0.03 (0.79) 
0.56 (0.07) 

0.86 

0.86 

0.35 


h -)■ 77/10"^ 
h -^ WW 

H ^ Tf 
A^Tf 

H^ -^ tb 


12.1 (2.1) 

0.36 (0.05) 

0.14 

0.14 

0.64 



Note that the benchmark point H has nonstandard 
Higgs bosons that are too hght to allow decays into hh. 
However, given that h is rather light in the region with 
suppressed bh couplings, it is possible that such exotic 
channels might be open, while still having an interesting 
di-photon signal. As mentioned before, in such cases it 
is possible that other channels involving SUSY particles 
are also open. 



CONCLUSIONS 



TABLE VIII. Masses and branching fractions in the BMSSM 
(and m the SM for h) for point H. 



tan/3 region of the MSSM. The charged Higgs can be 
searched for in the tUt channel. However, we emphasize 
again that the observation of light SUSY signals would 
give compelling evidence for BMSSM physics. There are 
also Tevatron covered (red) points at smaller mh values, 
around 110 GeV with an enhanced di-photon signal. We 
discuss these type of scenarios in the next section, to- 
gether with the Tevatron uncovered (blue points) in the 
same region. 



E. LHC searches 

Referring to Fig. [9] we split the Tevatron uncovered 
(blue) points according to whether their signal into pho- 
tons (and W's) is enhanced or suppressed. For the latter 
case, one has that h decays mainly into bottom and tau 
pairs. In these scenarios, h can be within the reach of 
the LHC in the tt channel, and if the suppression of the 
ZZ coupling is not extreme, maybe also in the ZZ — >■ Al 
channel. Higgs decay chains, such as i/ ^- /i/i —>■ 66rf, 
can also give rise to interesting (if challenging) signatures. 
We do not show a benchmark point here since the branch- 
ing fractions of the relevant Higgs decay chain modes will 
depend on the details of the sparticle spectrum. 

Point H: SM-like Higgs with enhanced di-photon signal 
We illustrate the features of models with a strong en- 
hancement of the di-photon signal with point H [? ] 
(shown in Table Villi. We see that h is rather light (the 
1/M^ operators contribute about 10 GeV to m^,), but 
escaped detection at LEP due to the strong suppression 
of the hh channel. The gg ^- h ^ 77 signal is larger than 
the SM one by a factor of 8, thus allowing for a very nice 
and clean detection of fi at the LHC. As was discussed in 
the context of Fig. [9] the same enhacement also occurs 
for the WW and ZZ channels. Therefore, and in spite of 
such a light Higgs mass, the gg -^ h ^ ZZ — > 41 channel 
would be at the reach of the LHC. 

For the remaining neutral Higgs bosons {H and A), 
one will have to consider the tt search. The charged 
Higgs may be detected at the LHC in the Th'r channel. 



We studied the Higgs collider phenomenology of 
BMSSM scenarios, i.e. supersymmetric extensions of the 
MSSM within an EFT framework where the effects of 
the BMSSM degrees of freedom enter through higher- 
dimension operators. As emphasized in [81 the first two 
orders in the 1/M expansion can be phenomenologically 
significant, and should be included. In the present work, 
we have performed a model-independent study to high- 
light the variety of collider signals that become available 
in such scenarios. 

The coupling of the lightest CP-even Higgs to bottom 
pairs can be suppressed due to cancellations between the 
MSSM contribution and those from the higher-dimension 
operators. It does not seem to require a special tuning 
of parameters and occurs in both the low and large tan /3 
regimes. As a result, the signals in clean channels, such 
as the di-photon or WW ones, can be greatly enhanced. 
This suppression in the tibb induces an enhancement in 
the gluon fusion production cross section, beyond the one 
arising from light sparticles in the loop. 

To emphasize the interplay between the Tevatron and 
the LHC, we have analyzed projections for the Tevatron 
assuming a total integrated luminosity of 10 fb~ per ex- 
periment and a 50% efficiency improvement in the WW 
and bb search channels with respect to present results. 
We find that the current Tevatron data already probes a 
large class of SUSY models, especially in the WW chan- 
nel. The future projections indicate that the bb channel 
can become effective for a SM-like Higgs search. More- 
over, a combination of the bb and WW search channels, 
together with the tt decay mode in the large tan/3 re- 
gion, would further enlarge the set of BMSSM models 
that can be probed at the Tevatron. However, our main 
interest in this work was to survey the types of signals 
that might be expected in SUSY scenarios, many of which 
are not realized in the MSSM limit. Improving the anal- 
ysis by combining channels and/or moderately increasing 
the luminosity will not significantly change our conclu- 
sions. Lightest CP-even Higgs bosons with masses above 
180 GeV, that can not be probed by the Tevatron, will 
be at the reach of the LHC. 

Most of the changes in the expected Higgs signals, com- 
pared to the MSSM, can be understood in large part 
from the altered Higgs spectrum. We have surveyed a 
wide range of possibilities by scanning over the parame- 
ter space of the higher-dimension operators. Motivated 
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by naturalness arguments, we have chosen the SUSY- 
breaking scale close to the EW scale, with the BMSSM 
physics at the TeV scale. In this case, the contributions 
from the SUSY particles to the Higgs spectrum are sub- 
leading compared to the ones coming from the BMSSM 
physics. In the case of a heavier SUSY spectrum, and 
for a scale M such that the effective field theory ap- 
proach remains valid, the qualitative features of the Higgs 
phenomenology triggered by the BMSSM physics will be 
similar. However, a detailed study should be performed 
for each specific choice of the heavy scale M, the scale 
of SUSY-breaking ms, and the /i-term to address the 
quantitative features of the Higgs sector. 

We have defined a number of "benchmark points" in 
order to discuss the correlations between different Higgs 
signals. Interestingly, we find that there can be signifi- 
cant mixing in the CP-even Higgs sector, allowing non- 
negligible couplings of both CP-even Higgs eigenstates to 
the EW gauge bosons. In addition, they can both be in 
the right mass range to decay predominantly into W's 
or Z's, thus enabling a detailed and direct study of the 
physics of EWSB. For these benchmark points, the 1/M^ 
effects add a few tens of GeV to mh, and have a rather rel- 
evant impact on the collider phenomenology (but we re- 
mind the reader that the 1/M^ operators can easily give 
a much larger contribution to 111^; see Fig. 1 of Ref. [5]). 
Furthermore, we have found viable examples where the 
nonstandard CP-odd Higgs can be produced in charged 
Higgs decays. Moreover, unusual decay chains such as 
h —>■ AA or _ff — >■ AA are also possible, without A being 
ultralight. These channels are most interesting in the low 
tan /? region where the tan /3-enhanced production of the 



nonstandard Higgs bosons is not available. These Higgs 
decay chains open the possibility of fully reconstructing 
the Higgs content of a 2HDM in such supersymmetric sce- 
narios. We also find scenarios where observing the Higgs 
sector is more challenging, and would require dedicated 
studies that go beyond the scope of this work. 

In conclusion, we find that Higgs signals in supersym- 
metric scenarios can be markedly different from those in 
the MSSM paradigm. If all third generation squarks turn 
out to be light {ms < 300 GeV), given the LEP Higgs 
mass bounds, this will imply a clear case for BMSSM 
physics. The heavier degrees of freedom could be at the 
kinematic reach of the LHC, but depending on their na- 
ture the direct discovery might be elusive. In either case, 
supersymmetric Higgs searches can provide evidence of 
physics beyond the MSSM. 
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